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ABSTRACT 

The galaxy cluster IE 0657-56 has a bullet-like subcl uster that is moving away 
from the centre of the main cluster at high speed. iMar kevitch et all l)2004j) 
recently estimated a relative velocity of Vbuiict = '^^OO^.Iqq'^ kms~^, based on 
observation s of the bow shock in front of the subcluster. The weak lensing 
analysis of IClowe et al.l l)2004j) indicates that a substantial secondary mass 
peak is associated with this subcluster. We estimate the likelihood of such a 
configuration by examining the distribution of subhalo velocities for clusters 
in the Millennium Run, a large ACDM cosmological simulation. We find that 
the most massive subhalo has a velocity as high as that of the bullet subcluster 
in only about 1 out of every 100 cluster-sized halos. This estimate is strongly 
dependent on the precise velocity adopted for the bullet. One of the ten most 
massive subhalos has such a high velocity about 40% of the time. We conclude 
that the velocity of the bullet subcluster is not exceptionally high for a cluster 
substructure, and can be accommodated within the currently favoured ACDM 
comogony. 

Key words: galaxies: clusters: individual (1E0657-56) - cosmology: dark 
matter - galaxies: formation 



1 INTRODUCTION 

Galaxy cluster IE 0657-56 is a singular example of a 
merging system. X-ray observations of the cluster have 
revealed a bow shock propagating in front of a bullet-like 
gas cloud moving away from the core of the main cluster. 
Based on the gas densit y jump across the shock front, 
IMarkevitch et al.l i2002r first derived an estimate of ve- 
locity of the bullet subcluster, Vbuiiet — 4000 km s~^. 
The uniq ue geometry of the system was recently ex- 
ploited bv lClowe et alJ (|200^, who combined the X-ray 
observations with weak lensing analysis to show that 
the mass distribution does not follow that of the hot 
cluster gas, but could plausibly be associated with a 
coUisionless dark matter component. These authors ar- 
gue that the data are not easily reconciled with mod- 
ified Newtonian Dynamics dMUgromlFigSal. MOND) in 
which the mass budget w ould be dominated by the gas. 
IMarkevitc h et al.l (|2Qq3) also use IE 0657-56 to derive 
constraints on the cross-section for self-interaction of 
the dark matter based on the spatial offset between the 
peaks in the X-ray and mass distributions. 

In this paper we estimate the probability of finding 
such a high velocity subcluster in a ACDM cosmology 
by examining the statistics of dark matter substructure 
halos (subhalos) in a very large cosmological N-body 
simulation. We calculate the distribution of subhalo ve- 



locities relative to their host halos for a large sample of 
cluster-sized hosts and use this to determine the fraction 
of clusters which contain a high velocity subcluster. 



2 PROPERTIES OF THE BULLET 
CLUSTER 

The cluster IE 0657-56 was discovered as an extended 
source in the Einstein imaging proportional counter 
(IPC) dat abase and was ident i fied as a rich cluste r of 
galaxies bv lTucker et alj lll995l) . lTucker et alJ lll99d) es- 
tablished the cluster redshift as z = 0.296 and identified 
the bullet subcomponent in a ROSAT X-ray image of 
the cluster. They also measured the temperature of the 
hot cluster gas as kT ~ 17 ke V, making IE 0657-56 one 
of the hottest known clusters. IMarkevitch et alJ i2002fl 
subsequently revised this to an average cluster temper- 
ature of 14-15 keV, albeit with large spatial variations, 
based on Chandra observations of the system. 

The transverse separation between the bullet sub- 
cluster and th e centre of the main cluster is ~ 
0.48 h-^Mpc iClowe et al.l |20o3). IMarkevitch et al.l 
||2004) estimate the velocity of the bullet subcluster, 
H>uiict = 4500lgjQ°kms~^; here and elsewhere we adopt 
a Hubble constant oi Ho = 73 kms~'^/Mpc. The line-of- 
sight velocity of the subcluster is ~ 600 kms~^ relative 
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to the cluster centre JBarrena et al.ll2002r) . so the direc- 
tion of the bullet's motion is very nearly in the plane of 
the sky. 

IClowe et al.1 (1200411 use weak lensing to model the 
mass distribution of the main cluster and the bullet sub- 
cl uster. They find that the main cluster is well fit by 
a iNavarro et al.1 l|l99fl . NFW) profile with concentra- 
tion C200 = 3.0 and virial radius r2oo = 1-64 /i~^Mpc, 
corresponding to a virial mass and velocity of M200 ~ 
2.16 X 10^5 /i-ijvfg and I/200 = (GM2oo/r2oo)'/' = 
2380 kms~^, respectively. They also find a secondary 
peak in the mass distribution that is clearly associated 
with the bullet subcluster, and they measure a mass of 
(5.3 ± 1.5) X 10^^ h'^Mo) for the subcluster within a 
cylinder of radius 0.11 /i~^Mpc. We adopt these values 
for the mass of the main cluster and the bullet rather 
than estimates based on the cluster galaxy velocity dis- 
persion or the gas temperature since the latter rely 
on assumptions of isotropy and hydrostatic equilibrium 
which may not be valid for a merging cluster system. 
We note that the virial mass we adopt corresponds to an 
overall X-ray temperature of kT ~ 15 — 1 6 keV accord- 
ing to the mass-temperature relations of lArnaud et al.l 
(|2005r . in good agreement with the observed X-ray tem- 
perature. 

To summarize, the main properties of this system 
relevant to this study are as follows: 

• Cluster IE 0657-56 is a very massive cluster, as 
evidenced by its high temperature and confirmed by the 
high mass detected by weak lensing analysis. 

• The bullet subcluster is the most massive substruc- 
ture in the cluster, and represents a mass ~ 2% that of 
the main cluster. 

• The bullet subcluster is moving away from the main 
cluster at velocity Vbuiict — 4500 kms~^ = 1.9 V200 and 
is at least 0.48 /i~^Mpc ~ 0.3 r2oo from the cluster 
centre. 

In the following section we assess the likelihood of 
such an object by examining the statistics of dark mat- 
ter halos in a large cosmological simulation. 



3 SEARCHING FOR THE BULLET 



This study makes 



use of the Millennium Run 
a very large cosmological N-body 

In 



arineel et alJl20 
simulation carried out by the Virgo Consortium.^ '^ 
this simulation a flat ACDM cosmology is adopted, with 
^dm = 0.205 and Q.h ~ 0.045 for the current densities in 
cold dark matter and baryons, h = 0.73 for the present 
dimensionless value of the Hubble constant, erg = 0.9 
for the rms linear mass fluctuation in a sphere of radius 
8 h~^Mpc extrapolated to 2: = 0, and n = 1 for the 
slope of the primordial fluctuation spectrum. The simu- 
lation follows 2160'^ dark matter particles from z = 127 
to 2; = within a cubic region 500 /i~^Mpc on a side. 
The individual particle mass is thus 8.6 x 10* h~^MQ, 



-*- http://www.virgo.dur.ac.uk/ 

^ http://www.mpa-garching.inpg.de/galform/virgo/millennium/ 



and the gravitational force is softened with a Plummer- 
equivalent comoving softening of 5 h~^kpc. Initial con- 
ditions were generated using the Boltzmann code CMB- 
FAST (Seljak and Zaldarriaga 1996) to generate a real- 
ization of the desired power spectrum which w as then 
impo sed on a glass-like uniform particle load iWhitd 
Il99el) . A modifie d version of the TREE-PM N-body 
code GADGET2 fSpringel et al.i'20011^,'2005') was used 
to carry out the simulation and full particle data are 
stored at 64 output times approximately equally spaced 
in the logarithm of the expansion factor. 

In each output of the simulation, halos are identi- 
fied using a friends-of-frie nds (FoF) groupfi nder with a 
linking length of 6 = 0.2 llDavis et alJll985ih . The virial 
radius, r2oo, is defined for each FoF halo by calculating 
the radius of a sphere, centered on the particle with the 
minimum potential, that encompasses a mean density 
200 times the critical value. Each FoF halo is decom- 
posed into a collection of locally overdense, self-bound 
substruc tures (or subhalos) us ing the SUBFIND algo- 
rithm of ISpringel et al.l (|2001a|). Of these subhalos, one 
is typically much larger than the others and contains 
most of the mass of the halo. We identify this as the 
main halo and subtract its centre of mass velocity from 
that of the remaining subhalos to compute the relative 
velocity between subhalos and their host halo. Here- 
after, we refer to this as the subhalo velocity, not to be 
confused with the internal circular velocity of the sub- 
halo. 

We search the Millennium simulation for halos 
whose most massive subhalo has a velocity relative to 
the main halo comparable to that of the bullet subclus- 
ter, i.e., l/sub > 1.9 1^200- We focus our search on cluster- 
and group-sized halos in the z — 0.28 output of the Mil- 
lennium simulation, the output closest in redshift to IE 
0657-56 . 

Figure shows the cumulative distribution of ve- 
locities for the most massive subhalo in host halos 
with M200 > 10" /i"^Mo, 3 X 10^" h-^MQ, and 
10^'' h~^MQ. The number of hosts in each of these 
mass ranges is A'hosts = 1491, 157, and 5, respectively. 
We note that IE 0657-56 is one of the hottest known 
clusters and that very few clusters in the Millennium 
Run have masses comparable to that of IE 0657-56 due 
to the limited volume of the simul ation. A ccording to 
the cluster temperature function of iHenrvl (|2p04), one 
expects les than 0.1 clusters as hot as IE 0657-56 in a 
volume similar to that of the Millennium simulation. In 
fact, at 2: = 0.28, the Millennium simulation contains 
one cluster halo with M200 > 2 x lO^'^ h~'^MQ. 

The shape of the velocity distribution appears rela- 
tively insensitive to host halo mass and the median sub- 
halo velocity is 1.1 ^200- We note that all subhalos have 
velocities have velocities much less than the maximum 
escape velocity, Hosc — 3.3 1^200 for an NFW potential 
with C200 — 6. 

The fraction of halos with A/200 > lO" H'^Mq 
whose most massive subhalo has a velocity greater than 
that of the bullet subcluster (scaled to the virial veloc- 
ity of the host halo) is 16 out of 1491 or approximately 
1%. However the velocity distribution drops steeply at 
high velocities, and this percentage increases (decreases) 
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Figure 1. Fraction of host halos whose most massive sub- 
halo has velocity greater than V^ub 7^200, for three different 
host halo mass ranges. The number of halos with M200 > 
10^* h-'^MQ, 3 X 10^* h-'^MQ, and lO^^ h'^MQ is 1491, 
157, and 5, respectively. The M200 > lO^* /i'^Mq distribu- 
tion is well fit by a eq. shown as the solid curve. The ver- 
tical lines indicate the velocity of the bullet subcluster (solid 
line) and lower and upper limits (dotted lines). About 1% 
of hosts have a most massive subhalo with velocity greater 
than that of the bullet. This percentage increases (decreases) 
to 10% (< 0.1%) for the lower (upper) limits of the bullet 
velocity. 



to 10% (-C 0.1%) if the lower (upper) limit is adopted 
for the bullet velocity. We no te that this is in a greement 
with the simulation results of lGill et alj ll2005l) who also 
find a small but significant fraction of high velocity sub- 
halos in cluster-sized host halos. 

We find that the velocity distribution is well fit by a 
function of the following form, shown as the solid curve 
in Figure m 



log 



iVi(>v;ub 

A^hosts 



V^sub/V200 
^10% 



(1) 



where i'io% is the velocity in units of V200 at which 
the fraction of halos drops to 0.1. Fitting this function 
to the M200 > 10" h'^MQ distribution at 2 = 0.28 
yields best fit values of 1110% = 1.55 and a — 3.3. We 
note that the best fit value of Wio% tends to decrease 
with decreasing redshift: Viq% — 1.64, 1.52, and 1.32 at 
z = 0.5, 0.11, and 0, respectively. The value of a shows 
no significant change with redshift and is typically con- 
sistent with a ~ 3.0 — 3.1. We attribute this to the 
increase in the virial velocity of the halo with respect to 
the velocities of subhalos at the time of infall. Indeed, 
the mean ^200 increases from 780 kms~^ at z = 0.5 to 
920 kms~^ at z = 0.0, whereas the mean subhalo veloc- 
ity decreases by only a few percent, from 1190 kms~^ 
to 1129 kms~^, over the same period. 

Having quantified the likelihood that a halo has 
a most massive subhalo with a velocity comparable to 



Figure 2. Fraction of host halos whose i-th most massive 
subhalo has velocity greater than V^ub/V'200, for 1st, 2nd, 
3rd and 10th most massive subhalos. More massive subhalos 
are biased toward lower velocities. Solid curves show fits with 
eq.0 with values for «io% of 1.55, 1.64, 1.72, and 1.79 in order 
of increasing subhalo rank, respectively. 



that of the bullet subcluster, we now estimate the proba- 
bility of finding a high velocity subhalo amongst a halo's 
n most massive subhalos. The probability of drawing at 
least one subhalo with velocity > 14ub from the n most 
massive subhalos is given by 



p{> v;,b) = 1 - n 



1- 



iV,(> 14ub) 

A'^hosts 



(2) 



We investigate whether subhalo velocity is indepen- 
dent of subhalo mass in Figure |21 by comparing the ve- 
locity distributions for the 1st, 2nd, 3rd and 10th most 
massive subhalos in halos of mass M200 > 10^* /i~^Mq. 
We find that more massive subhalos are slightly biased 
toward lower velocities. The best fit values of iiio% are 
1.55, 1.64, 1.71, and 1.79 for the velocity distributions of 
the 1st, 2nd, 3rd and 10th most massive subhalos, and 
i;io% — 1.8 for subhalos of higher rank. We find that the 
following formula accurately describes the trend: 



Vio%{i) = 1-8 - 0.25 exp(-0.45 (i - 1)). 



(3) 



Note that the value of a decreases slightly with in- 
creasing subhalo rank, but we find that a constant value 
of Q — 3.1 provides an adequate fit to the distributions 
for all subhalo ranks. 

We combine eqs.0|5| andj^in order to predict the 
probability of at least one subhalo with Vaub > 1.9 V2oo- 
This gives a probability of 39.7% which agrees well with 
the actual fraction, 40.8%, found for halos in the Millen- 
nium simulation with M200 > 10^"* h~^MQ. In compar- 
ison, if we adopt a constant value for Uio% of 1.55 (1.8) 
we predict a probability of 12.4% (49.3%). 

In order to convert the subhalo mass rank into a 
fraction of the host halo mass, we compute the sub- 
halo mass function for our halo sample. We find that 
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the differential subhalo mass function for halos with 
M200 > 10^'' h~^ Mq is well fit by a power law with 
slope —0.88, in agreemen t with the ma ss functions of 
cluster halos presented bv IPe Lucia et al. ( 2004) . Inte- 
grating the differential mass function yields the total 
number of subhalos with mass greater than Mg^b '■ 



logAf(>Af,ub) = -0.88 log 



Ms, 



M200 



1.67. 



(4) 



Solving this equation for Afsub gives the subhalo 
mass corresponding to subhalo mass rank A*'. For exam- 
ple, for A*' = 10 we find Afsub = 0.001 M200, correspond- 
ing to subhalos of mass lO^'^ h~^MQ for cluster-sized 
host halos of mass M200 — 10^^ /i~^A/q. Our previ- 
ous result therefore implies that 40% of massive cluster 
halos have at least one subhalo of mass greater than 
10^^ h~^MQ with a velocity comparable to that of the 
bullet subcluster. 

We now return to the bias in the velocity distri- 
bution of the most massive subhalos (see Figure 0. 
This is related to the fact that more massive substruc- 
tures are preferentially located in the external regions 

I ■ 1 I 1. 

of their host halos, as noted bv lDe Lucia et alJ Il2004li . 

These subhalos are closer to the apocentre of their or- 
bits, and therefore have lower velocities compared to 
subhalos near pericentre. Conversely, we expect to find 
high velocity subhalos near the centre of their host ha- 
los. Indeed, of the 16 high velocity subhalos that are 
bullet subcluster candidates, all but one are located at 
r < 0.6 T-200, whereas only 15% of total sample of most 
massive subhalos are found within this radius. This cor- 
relation between subhal o velocity and clu stercentric dis- 
tance was also noted bv lGill et alJ (|2003). 

In Figure|21we compare the velocity distributions of 
subhalos in the inner and outer regions. Subhalos within 
the central 0.6 7-200 of the host halo (about half of the 
total within r2oo) are indeed biased to higher velocities. 
The best fit values of Uio% are 1.74, 2.16, and 1.87 for 
the velocity distributions of subhalos with r < 0.6 r2oo, 
r > 0.6 r2oo and the combined sample, respectively. We 
also note that velocity distribution of outer subhalos is 
not as well fit by eq.^as the other distributions, however 
the deviations are typically < 0.1 dex. 

Finally, we note that the direction of the bullet 
subcluster's velocity is an additional constraint that we 
have not so far considered. As stated in ^ the bullet 
subcluster is moving away from the cluster centre, i.e., 
has a positive radial velocity. Of the 16 bullet subcluster 
candidates, five have positive radial velocities. Overall, 
we find that a smaller fraction of massive subhalos tend 
to have positive radial velocities compared to less mas- 
sive subhalos; this fraction is 20%, 30% and 40% for 
the most massive, second most massive, and tenth most 
massive subhalo samples. 

We attribute this to the mass loss and disruption 
of subhalos by tidal stripping, which occurs near orbital 
pericentre. This can d ecrease the mass of a subhalo by 
as much as 90% (e.g.. lllavashi et al.ll200a) . potentially 
downgrading the mass rank of more massive subhalos 
after they pass through pericentre. Less massive sub- 
halos can be disrupted altogether resulting in the de- 
pletion of subhalos with positive radial velocities. How- 
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Figure 3. Fraction of subhalos with velocity greater than 
V^ub/V200, at different radii. Subhalos in the outer regions 
are biased toward lower velocities. Solid curves show fits with 
eq.0 with values for 1110% of 1.74, 2.16 and 1.87 for subhalos 
with r > 0.6 r200i with r < 0.6 r2oo arid for the total sample, 
respectively. 



ever, the velocity distributions of subhalos with positive 
and negative radial velocities are very similar; subhalos 
with positive radial velocities are slightly biased toward 
lower velocities, but the difference in the best fit values 
of Wio% and a are typically ~ 0.05 and ~ 0.5, respec- 
tively. We therefore conclude that if the direction of the 
bullet subcluster is considered, the probabilities we have 
estimated are reduced by about 60 — 70%, the fraction 
of subhalos with negative radial velocities. 



4 CONCLUSIONS 

The bullet subcluster is a massive substructure in galaxy 
cluster IE 0657-56 moving at a high velocity relative to 
the centre of the main cluster. In units of the virial ve- 
locity of the cluster, the velocity of the bullet subcluster 
is Vbuiict — 1.9 V200. We have examined the distribution 
of subhalo velocities relative to their host halos in a large 
cosmological simulation in order to assess the likelihood 
of a system like IE 0657-56. 

We calculate the velocity distribution of the most 
massive subhalos in 1491 host halos with virial masses 
M200 > 10" H'^Mq and find that about 1 in 100 have 
velocities comparab le to that of the bullet s ubcluster if 
the best estimate of lMarkevitch et alJ l|2004|) is adopted 
for the bullet velocity. However, this fraction depends 
strongly on the velocity of the bullet and ranges from 
<^ 0.1% to 10% for the upper and lower limits on the 
bullet velocity, respectively. We find that more massive 
subhalos are biased towards lower velocities, as are sub- 
halos in the outer regions of halos. Taking this into ac- 
count, we find that at least one of the ten most massive 
subhalos has a velocity as high as that of the bullet sub- 
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cluster in 40% of all host halos. We also find that sub- 
halos are preferentially found to be moving toward the 
centres of halos, most likely a result of tidal depletion of 
subhalos at pericentre. With this additional constraint, 
the likelihood of the bullet cluster drops to about 1 in 
500. We conclude that the best estimate for the veloc- 
ity of the bullet subcluster is high but not extraordinary 
considering the mass of its host cluster. It is a rare but 
not an impossible event within the currently favoured 
ACDM comogony. 
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